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The properties of quark and gluon jets depend on jet definitions 
and event selection. I discuss how these can be included in cal- 
culations and present jet definitions designed to give unbiased 
jets. 



1 Introduction 



Perturbative QCD predicts the ratio of multiplicities in quark and gluon jets 
to reach an asymptotic value Ca/CV = 9/4. At accessible energies, it is 
however found to be substantially lower, ~ 1.5 at 90GeV.|l[ Instead of a 
method to establish the colour factor ratio, the average multiplicities are a 
simple and powerful probe of subleading corrections. Energy conservation 
effects on the multiplicity ratio have been studied|2|, fj] and results are in 
agreement with data. Q 

As there is no colourless point source for gg events available at energies 
above onium resonances, the major way to study scale evolution of N gg is 
via jets. However, event selection and jet definitions introduce a bias that 
needs to be well understood. That is the issue of this talk, with focus on 
e + e~ annihilation. 

2 Two-jet Events 

We start to investigate the bias in the simple case of two-jet events selected 
with a jet clustering algorithm with a resolution scale k± r . The situation is 
illustrated in Fig. ([]]), in the plane of logarithmic variables lnf^/A 2 ) and 
rapidity y. This coordinate system is well suited to illustrate the multiplicity 
formulae below, as the emission density of gluons is essentially flat in rapidity. 

In a central rapidity range the upper constraint on sub jet transverse mo- 
menta is k± r , but in a forward cone of each jet the kinematical constraint is 
more restrictive. Summing up the multiplicities in the forward cones and the 
central region, we getR ^[ 

, C = ei (1) 

Q 2 =C,k\ r 

The scale shift factor Q is a modified leading log (MLLA) correction to the 
leading order result ( = 

Fig. (0) shows the ratio biased over unbiased multiplicities for different 
resolution scales k^ r at fixed energy y/s = 90GeV. To calculate the biased 
multiplicity in Eq. ([[]) we have used a simple fit to Ariadne+Jetset MC0 
results for the unbiased A q q. As seen, the effect of the MLLA correction 
factor ( is significant. The bias is in general important, being more than a 
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Figure 1: In a two-jet event sample selected 
with a resolution scale k± r , gluon momenta 
are constrained by k±_ < k± r , in addition 
to the kinematical limit \k\ < \fsj1 which 
corresponds approximately to the triangle 
\y\ < ^\ia(s/kj_) . The multiplicity in the 
forward cones of the jets (dark gray) is un- 
biased, while the multiplicity in the central 
region (light gray) is constrained by k± r . 
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Figure 2: The ratio biased 
over unbiased multiplicity. 
The bias gives a correction 
larger than 10% for k± r < 
20GeV. The MLLA pre- 
diction is in agreement with 
MC simulations using the 
Durham jet finder. 



10% effect for k ±r < 20GeV. Also shown in Fig. 
Ariadne and the Durham jet algorithm. 



are MC results with 



3 Three-jet Events 

The multiplicity in gluon jets can be extracted from qqg three-jet events. 
If a three-jet event sample is selected with a fixed resolution k± r , we get 
complicated scale dependences, and all jets are biased. It is therefore more 
suitable to perform iterative clustering until exactly three jets remain, and 
study the event as a function of the k^ g of the softest jet, presumably the 
gluon jet. In this approach the gluon jet is essentially unbiased, but the bias 
in the q and q jets needs still to be considered. 

The qqg configuration is illustrated in Fig. @. The emission density of 
softer gluons corresponds to radiation from a qg- and a gq colour dipole.|| 
There is also a colour correction term, corresponding to a qq dipole, whose 
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Figure 3: The qqg configuration, (a) The emission density corresponds to 
two dipoles plus a colour suppressed correction here represented by a dipole 
spanned between the q and q. (b) The correction term can be approximated 
by a colour factor reduction N c /2 — ► Cp in the quark ends of the two dipoles 
attached to the gluon. (c) In the logarithmic phase space variables, the two 
dipoles can be represented by the original qq triangle plus a double sided 
fold. The magnitude of the Cp regions in the q and q jets introduces a new 
scale, here called E qq . 

contribution is weighted with the negative factor —N~ 2 . As apparent from 
the negative weight, this colour correction dipole is not an independent emit- 
ter. Instead it can be assumed to reduce the colour factor from N c /2 to Cp 
in the quark ends of the two other dipoles. The magnitude of the Cp phase 
space regions can not be determined by perturbative QCD, since the notion 
of an infrared cut-off independently defined for each dipole does not apply to 
the colour correction term. This introduces an uncertainty of relative order 
7^2 jj- to the total hadronic multiplicity and hinders the possibility to fully 
determine corrections of relative order 1/ In s from first principles. 

The magnitude of the Cp regions in the q and q jets introduces a new 
scale, here called £ qq . Provided E qq is not much smaller than s, the total 
multiplicity for the three-jet event is|| 

iV qqg « iV qq (£ qq , kl g ) + ± Ngg (k 2 *). (2) 

Z ^qq 

This equation shows the possibility to extract unbiased N gg using three-jet 
event data. The importance of taking the bias on the qq system into account 
is illustrated in Fig. (^). The result is sensitive also to the assumed value of £ 
in Eq. ([I]). Neglecting the MLLA corrections to the bias, by setting £ = 1 in 
Eq. ([[]), raises the prediction on iVg g by about two to three charged particles. 
This emphasizes the importance of confronting Eq. (|TJ) with data. 
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Fi gure 4: Ngg extracted from three-jet events with y/s — 90GeV and wSq^ — 
60GeV. Neglecting the bias (dashed line) underestimates N gg at low scales, 
as compared to the result when the bias is taken into account (solid line). 



Though the dependence on E qq - is formally suppressed by l/A^ 2 , different 
E q g values implicitly correspond also to different recoil assumptions. The 
recoil treatment when qqg — > qqgg can in principle be constrained by com- 
parisons with the four parton matrix element, [0] but an exact solution is 
beyond reach within a parton cascade formalism. This puts a limit on the 
accuracy in the analysis presented here. The prediction on N q ^ s for the rea- 
sonable assumptions £ q q = s qq - and £ qq - = s differ by about | charged particle 
when yfi = 90GeV and = 60GeV.[| 



4 Individual Jets 

The method described in the preceding section has the advantage that we 
need no explicit investigation of individual jets, but it may suffer from depen- 
dences on recoil assumptions. In this section we discuss the scale dependences 
in jets, and present a jet definition designed to give unbiased gluon jets. This 
gives an alternative way to observe the scale evolution of the unbiased iVg g . 

The properties of a jet depend on (at least) two scales. The maximum k± 
for subjets and the available rapidity range Yj (at some cut-off scale A) in the 
jet, which depends on the jet energy. The two scales are truly independent, 
and can not be combined into one single effective scale. As an example, each 
hemisphere of a two-jet event illustrated in Fig. (|I|) can be seen as a quark 
jet with transverse momentum scale k± r and a rapidity scale |ln(s/A 2 ). 

In general, Yj is not easily determined. Each jet can get multiplicity con- 
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Figure 5: The "Boost" jet 
boundaries. The gluon jet 
is unbiased and the quark 
jets have well defined rapid- 
ity scales Yq. 



Figure 6: MC simulations with the Boost 
method. For fixed k± g the charged multi- 
plicity in the gluon jet is essentially inde- 
pendent of y and in good agreement with 
MC simulations of \N gg . 



tributions from several dipoles, these dipoles are boosted away from their rest 
frame, and the opening angle of the jet need not be azimuthally symmetric. 
In order to get better control of the Yj scale, we suggest a jet analysis in two 
steps 

• Find jet directions using a fc^-based cluster algorithm. 

• Redefine jet boundaries to get well controlled Yj values. 

With this approach it is possible to define jet regions as illustrated in Fig. 
This implies that we can study unbiased gluon jets and also quark jets with 
well defined rapidity scales 

K q(g) = iln( S ) +(-)y, (3) 

where y = | ln[(l — Xq)/(1 — x q )} is the rapidity of the gluon jet. 

The boundaries are defined in terms of a Lorentz boost which is described 
in detail elsewhere. Fig. ([]) shows the results on N g from MC simulations 
using this "Boost" method. The multiplicity is plotted as a function of k± g , 
but for each value of k± g results for several values of E g are shown. The 
crosses are mostly on top of each other, confirming that the method gives 
unbiased gluon jets depending on one scale only. 

The Boost method enables also a study of the F q dependence of the biased 
multiplicity for fixed k± g , complementary to the study suggested in section |2|, 
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Figure 7: MC simulation of quark jets with the Boost definition (crosses). 
For fixed k± g , the y dependence is linear and in agreement with Eq. (|]) (solid 
line) . 

where the k±_ r dependence for fixed "rapidity scale" In s is studied. Fig. (|7|) 
shows that MC results for the quark jets in the Boost method agree with the 
analytical expression in Eq. ([]]). 

5 Summary 

Event selection introduces bias on multiplicities in two-jet and three-jet 
events. An expression for this bias,|| valid in the modified leading log ap- 
proximation, awaits confrontation with data. Taking the bias into account, 
the unbiased multiplicity iV g g can be extracted as 2(iV qq - g — iV qq -(s, &i g )).fj 

In general, jets are biased, and their properties depend on two scales, 
transverse momentum and an available rapidity range. Unbiased gluon jets 
where the two scales coincide can, however, be defined by the "Boost" bound- 
ary definitions. || The Boost method also gives quark jets where both scales, 
though different, are well defined. An expected linear dependence on the 
rapidity scale is seen in MC simulations. 
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